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The MINOS experiment at Fermilab has recently reported a tension between the oscillation 
results for neutrinos and anti-neutrinos. We show that this tension, if it persists, can be 
understood in the framework of non-standard neutrino interactions (NSI). While neutral 
current NSI (non-standard matter effects) are disfavored by atmospheric neutrinos, a new 
charged current coupling between tau neutrinos and nucleons can fit the MINOS data without 
s ! violating other constraints. In particular, we show that loop-level contributions to flavor- 

(3JT) violating r decays are sufficiently suppressed. However, conflicts with existing bounds could 

arise once the effective theory considered here is embedded into a complete renormalizable 
model. We predict the future sensitivity of the T2K and NO^A experiments to the NSI 
parameter region favored by the MINOS fit, and show that both experiments are excellent 
tools to test the NSI interpretation of the MINOS data. 

i 1 PACS numbers: 14.60.Pq, 12.60.-i, 13.15.+g 

Oh 

1. INTRODUCTION 

i—i Recently the Fermilab MINOS experiment has reported new results on £v disappearance [lj. 

Interestingly, the values of the neutrino oscillation parameters Am| 2 and sin 2 2#23 preferred by 
this anti-neutrino measurement are in tension, at the 90% confidence level, with the preferred 
region in the Am 2 2 versus sin 2 2023 plane for the neutrino disappearance. A likely explanation 
of this tension is lack of statistics especially in the anti-neutrino channel where the number of 
observed events in the far detector is approximately 100 whereas for the neutrino channel the 
number of events in the far detector is of order 20 times larger. One could speculate that the 
tension between neutrino and anti-neutrino disappearance might be the first hint of CPT violation 
in the neutrino sector pHS]. However, given that CPT conservation is such an important tenet 
of modern quantum field theory, it is important to explore other possibilities for new physics 
before giving up CPT conservation. Since the MINOS experiment is performed not in vacuum but 
with more than 700 kilometers of Earth matter between the source and the detector, there is the 
possibility of Wolfenstein type matter effects [B] leading to apparent CPT violation. While matter 
effects are not relevant to za, disappearance in MINOS in the standard three-flavor framework, they 
may become important in scenarios with sterile neutrinos [7], or if new non-standard interactions 
are contributing to the potential that the neutrinos experience when traveling through matter. 
The latter option will be explored in this paper. In addition, we will also consider non-standard 
interactions modifying the neutrino detection process in a CP non-conserving way. 

Non-Standard Interactions (NSI) in the neutrino sector have been introduced first as an al- 
ternative to standard oscillations [6], and later as a possible addition IS]. The phenomenology of 
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such subdominant NSI effects in current and near-future accelerator neutrino experiments has been 
investigated by many authors [51-I16|. and bounds on NSI have been derived from oscillation and 
non-oscillation data |17H20j . In the context of the latest MINOS results, NSI have been brought 
up in [2H [22] , and a concrete model has been proposed in [23] . 

In the following, we will first show analytically how NSI affect neutrino oscillations in the two 
flavor limit (sec. [2]), and then perform fits to the MINOS and data including different types of 
NSI (sec. [3]). In sec.|4j we will consider the potential of the T2K [21] and NO^A |25j experiments to 
test the NSI interpretation of MINOS. Finally, we will discuss our results from the model-building 
point of view and draw our conclusions in sec. [5] 



2. ANALYTICAL FRAMEWORK OF NON-STANDARD INTERACTIONS 



2.1. Neutral current NSI 



At O(GeV) energies relevant to neutrino oscillation experiments, non-standard interactions can 
be introduced in the Lagrangian as effective dimension 6 operators coupling neutrinos and charged 
fermions. We will first discuss new neutral current couplings of \x and r neutrinos to normal 
matter, i.e. electrons, up-quarks, and down-quarks. Phenomenologically, such operators will result 
in non-standard matter effects, so that the Hamiltonian governing neutrino propagation in the \i—t 
sector will read 



\ 1 f m f m 



(i) 



where 



_ / cos6>23 sin 6*23 \ ^ 

\- Sin 6*23 cos 023/ 

is the leptonic mixing matrix with the mixing angle 023, E is the neutrino energy, and A = 
2\/2GfN £ E is the matter potential depending on the electron number density N e along the neu- 
trino trajectory. The parameters ej^, , and give the relative strength of the non-standard 
interactions compared to Standard Model weak interactions. The superscript m indicates that 
these parameters describe non-standard neutrino matter effects. e™ T can in general be complex, 
while ej^ and have to be real in order to preserve the hermiticity of the Hamiltonian. In 
the following, we will set e™^ = since terms proportional to the identity matrix do not affect 
the neutrino oscillation probability, implying that oscillation experiments are only sensitive to the 
combination e™ — e^. 

The disappearance probability for v„ in matter of constant density is given by 

/Am 2 L\ 

P(vp -+ = 1 - sin 2 20jv sin 2 f j . (3) 

where 0jv and Am 2 N are the mixing angle and mass squared difference in matter. These matter 
oscillation parameters can be found by solving the following set of coupled equations, 

Am 2 v cos20 Ar = Am| 2 cos 20 2 3 + e™ T A , (4) 
Am% sin 207V e i<t>N = Aml 2 sin20 23 + 2e™ A . (5) 
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The phase (f>N is unobservable in oscillations and can be phased away. The solution is trivial to 
find and is given by 



Am 2 N = ^(Am 2 32 cos26 2 3 + e™A) 2 + \Am 2 32 sm2e 2 3 + 2e™A\ 2 (6) 

and sin 2 2<9jv = I Am§ 2 sin 26 23 + 2e™ A\ 2 / Am% . (7) 
Thus, the disappearance probability can be written as 



|Am^ 2 sin2g 23 + 26^| 2 . 2 ( Am 2 N L \ 
P(v^u») = 1 ^ sm 



(8) 



with Am 2 N given by eq. ([6]). This expression agrees with that found in ref. p3] when expanded to 



first order in e™ T and e™ T . 



In vacuum, ^4 = 0, and eq. (|8|) reduces to the standard two flavor disappearance probability. In 

2 N L/AE) « (Am 2 N i 

2e™A 



the small L/E limit, i.e. when sin 2 (Am 2 v L/4£') « (Am%L/4E) 2 , 



sin 26>23 + 



Am 2 2 



Am§ 2 £ 
4£ 



(9) 



so ej^. modifies the disappearance probability in this limit whereas does not. 
For anti-neutrinos, e™ T — > e™* and A — > —A, so that in matter 

-> ^) ^ -> *m) (10) 

without CPT violation. 

We note three interesting symmetries in eq. Q: First, the expression depends only on 
cos[arg(e™ )], not on sin[arg(e™)]. Therefore, it is invariant under the replacement 

arg(e™)^27m-arg(e™) (11) 

for arbitrary integer n. Moreover, it is easy to verify that Piy^ —> v^) is also invariant under the 
simultaneous replacements 

e^-e™, Am 2 2 ^-Am 2 2 . (12) 



as well as under the transformation 



#23 -> \ ~ 023 ■ (13) 



These symmetries will generate an eightfold degeneracy. For fixed E, P(vn v a) has an additional, 
continuous, symmetry: It is invariant under any simultaneous variation of Am 2 2 , #23 , > l e J^-|> 
and arg(e^) that leaves Am 2 N and sin 2 2#7v invariant. If we demand this invariance for neutrinos 
and anti-neutrinos, we obtain 4 equations for 5 free parameters, implying that the symmetry is 
continuous. However, since A is energy-dependent, this symmetry will not be manifest in the 
MINOS data, which covers a broad range of energies, and we will therefore not consider it further 
in this paper. 

Note that these symmetries are exact only in the two-flavor framework. In the three-flavor case 
with large #13, it is for example possible to determine the mass hierarchy by observing standard 
matter effects either in the z/„ — > v e channel or directly in the v„ — > channel |26j . 

In fig. [T] we have plotted the disappearance probabilities with all three combinations of non-zero 
and/or e™ for representative values of the e m 's for the MINOS experiment. One can see that 
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Figure 1: The survival probability for MINOS (L — 735 km) with a variety of neutral current NSI (non- 
standard matter effects) turned on as indicated in the plot. The solid (blue) lines are the neutrino survival 
probabilities whereas the dashed (red) lines are for anti-neutrinos. The dotted (black) lines are the vacuum 
survival probabilities. For the standard oscillation parameters, we have assumed Am§ 2 = +2.86 x 10~ 3 eV 2 
and sin 2 (9 23 = 0.38. 

non-zero e™ r changes the disappearance probability most notably at large energies and shifts the 
position of the minimum in energy. Whereas non-zero changes the disappearance probability 
most notably near the first oscillation minimum, especially in the depth of the minimum. Since the 
tension between MINOS neutrino and anti-neutrino data is both in the position of the minimum 
and in its depth, one requires non-zero e™ r and non-zero e™ in order to lift the tension in the 
optimal way. 

2.2. Charged current NSI 

As an alternative to neutral current NSI, we also discuss non-standard charged current inter- 
actions affecting the neutrino production and/or detection processes as an explanation for the 
MINOS results. If the Wilson coefficients of the corresponding effective operators are complex, the 
interference term between the standard and non-standard Feynman amplitudes can be different 
for neutrinos and anti- neutrinos and CP-violating phenomena can emerge. The modifications to 
the far detector event spectra observed in MINOS can be induced by (i) operators leading to a 
modified flux of at the far detector, but not at the near detector, and (ii) by operators leading 
to the production of muons in interactions of v T . (We neglect the possibility of non-standard inter- 
actions of v e since their flux at the far detector is between one and two order of magnitude smaller 
than that of Vn because of the low v e contamination of the NuMI beam and the smallness of the 
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mixing angle #13.) The only way of realizing case (i) in a three-flavor framework is to postulate 
a v T contamination in the NuMI beam, which would be invisible to the near detector, but would 
have partly oscillated into f« when reaching the far detector. However, results from the NOMAD 
experiment [27j constrain the v T contamination of a NuMI-like neutrino beam at short baseline 
to be less than 1.7 x 10~ 4 at 90% confidence level, much too small to be relevant to MINOS. We 
therefore neglect case (i) in the following, and focus on case (ii), a non-standard interaction of the 
type 

v T + N -> X + n , (14) 

where N is a nucleon and X stands for the hadronic interaction products. The operator generating 
this process has the structure 

£nsi 3 -2V2G F €%V ud [uYd] \flj p P L u T ] + h.c. , (15) 

where u and d denote the up- and down-quark fields, Gf is the Fermi constant, Pj, = (1 — 7 5 )/2, 
and gives the strength of the non-standard interaction compared to Standard Model weak 
interactions. (In the notation from ref. [20J, this coefficient would be called .) In principle, 
one could also consider the axial- vector operator [u^ p ^d] [p, r ypPLi' T ], but, being parity-odd, this 
operator would lead to the decay ir — > ^lv t1 which is strongly constrained by NOMAD, as discussed 
above. We will therefore neglect axial-vector NSI. We will also not consider scalar, pseudo-scalar, 
and tensor operators, since they could lead to the required interference between standard- and 
non-standard Feynman amplitudes only if the outgoing muon in the detector undergoes a helicity 
flip [14J. For O(GeV) muons, this would correspond to a suppression of the non-standard interaction 
rate by mJ 1 E ~ 0.1. 

If is non-zero, the counting rate in MINOS is no longer proportional to the standard survival 
probability P{yu, v u)-, but rather to an apparent survival probability P(vn — > Vn)-, defined by 
the number of muons produced in the detector in interactions of neutrinos with a given energy E, 
divided by the number of muons that would be produced in the absence of neutrino oscillations and 
non-standard interactions. Thus, P — > Vn) includes the possibility that the neutrino flavor has 
changed into v T during propagation, but non-standard interactions still lead to a final state muon, 
normally associated with interactions. Since the amplitudes for the processes + N — > X + fi 
and Vp 2^,. Vr _|_ _/y _^ x + /i can interfere, P{y^ — > Vp) can be larger than unity and is therefore 
not a survival probability in the usual sense. 

In the two-flavor approximation we find for -P(f^ —5- 



|2 



1 + 2 |e*J cot 20 2 3 cos [ arg(e^)] - \e d T 

Am?> L\ ( AmloL 



sin 2 2023 sin 2 | — r3 



AE 



P(^ -> I'm) = 1 

+ 2|4|sin2023 S m[arg(4)]sin( A ^)cos( A ^| L ) . (16) 

For anti-neutrinos, the sign of arg(e^) has to be reversed, and thus 

P(^ -> ^) + Pfa -> ^) • (17) 

without CPT violation. 



Note that the non-standard terms proportional to |e^| and (e^l 2 in the first line of eq. ( 16 ) have 



the same energy dependence as the standard oscillation term and can therefore change only the 
depth of the oscillation dip, but not its position. They can thus only change the apparent value of 
sin 2 2023 that would be reconstructed in a standard oscillation analysis neglecting NSI. Moreover, 
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Figure 2: The apparent survival probability for MINOS (L = 735 km) without NSI (black dotted curve), and 
in a scenario with non-zero e^ M . The solid (blue) line in this case is for neutrinos, while the red (dashed) line 
is for anti-neutrinos. For the standard oscillation parameters, we have assumed Am| 2 = +2.74 x 10 -3 eV 2 
and sin 2 9 2 3 = 0.41. 



these terms do not depend on the sign of arg(e^) and can therefore not introduce an asymmetry 
between neutrinos and anti-neutrinos. The interference term between standard and non-standard 



amplitudes in the second line of eq. (16), on the other hand, can be different for neutrinos and anti- 



neutrinos. Since eq. (16) is invariant under the simultaneous replacements arg(e^„) — > — avg(ef^) 



and Ain^ 2 L/AE — > ir — Am^L/AE, the depth of the oscillation minimum will be the same for 
neutrinos and anti- neutrinos, and only its position will be different. This is also illustrated in 
fig. 2j where we plot Piy^ — > v^) including non-zero ef^ for neutrinos and anti- neutrinos. 



ike eq. <\8ti for neutral current NSI, also eq. (16) exhibits several symmetries. In particular, 



the expression is invariant under the simultaneous replacements 



arg(e T 



27rn — arg(e 



d ) 



Am§ 2 



-Am| 2 



(18) 



and under the transformation 



arg(el) -> (2n + l)ir - arg(e^. 



#23 



7T 

2 



#23 



(19) 



for arbitrary integer n. Actually, the second of these symmetries can be generalized to a contin- 



uous symmetry. To see this, note that eq. (16) is invariant under simultaneous changes of \e d 



Tfli 



arg(e^), and 623, provided that the coefficients of the energy dependent factors sin [Am\ 2 L / AE] 
and sm[Am^ 2 L / AE] cos[Am^ 2 L / AE] remain invariant. This requirement imposes two constraints 
on the three parameters |e!^|, arg(e^ M ), and 6*23, so that there will be an infinite set of solutions. 
Note that the symmetries of eq. ( |16[ ), like those of eq. ([8]), are exact only in the two- flavor frame- 
work. 
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3. FIT TO MINOS DATA 

To test the compatibility of MINOS data with the hypothesis of non-standard neutrino inter- 
actions and to extract the allowed values for the NSI parameters, we have performed fits using 
a modified version of GLoBES [28, 29J (see appendix [A| for details). We have first checked that 
we were able to reproduce the results of the two-flavor standard oscillation fits performed by the 
MINOS collaboration. We find that the fit to only (z^) data yields x 2 / dof = 12.3/12 (2.3/3), 
while a combined fit to both data sets results in x 2 /dof = 20.1/17. According to the param- 
eter goodness-of-fit test described in ref. [30], this means that the probability for the apparent 
inconsistency between the and data sets to be merely a statistical fluctuation is about 6%. 
Even though this probability is still relatively large, let us now study how the fit improves if the 
possibility of non-standard interactions are included. 



3.1. Neutral current NSI 



We begin by assuming only neutral current NSI (see sec. 2.1 ). Using a three-flavor fit including 
the NSI parameters |e™ |, arg(e™ ), and , we find the following parameter values for the best fit 
point: 



e™ = -0.40 = 0.40 e 1,0Mr sin 2 # 23 = 0.38 

e™ = -2.16 Am 2 32 = +2.86 x 10~ 3 eV 2 ^ 



The x 2 a t the best fit point is 12.6, and the number of degrees of freedom is 14. Note that the 
best fit value for e™ T is essentially real and negative. An equally good fit is obtained when the sign 
of arg(e^) is inverted (eq. (11)), when the signs of e™ T , e™, and Am| are flipped simultaneously 



(eq. (12)), or when the sign of e™ and the octant of 623 are changed (eq. (13)). This means 
that three-flavor effects are not large enough to spoil the symmetries of the two-flavor survival 
probability eq. Q. Thus, the best-fit point is eightfold degenerate. In fig. [3j we compare the 



theoretically predicted event spectrum at the parameter point eq. (20) (blue dotted histograms) 
to the MINOS data and to the spectra obtained from a two-flavor standard oscillation fit. 

Our fit is in qualitative agreement with the one in ref. |21j . but points to somewhat larger 
values of |e^_| ~ 0.4 compared to |e™ | ~ 0.1 in ref. [21]. We surmise that the main reason for this 
disagreement is that our fit is based on a simulation of the MINOS event spectrum, while, to our 
understanding, the authors of ref. |21| have directly fitted the oscillation probability to the ratio of 
observed to expected event numbers, as computed by the MINOS collaboration. A fit at the level 
of the oscillation probability, however, cannot fully include experimental energy resolution effects. 
The fit to the spectrum used in our work is able to reproduce the MINOS best fit points as well 
as the allowed regions in the Am^-sin 2 2023 plane rather well (see fig. [8] in appendix |a|). 

The allowed regions and best fit values for |e^|, arg(e^j-), and |e™| are shown in the three 
panels of fig. |4j We find that, in spite of the large best fit value for |e™ |, MINOS is compatible 
with = at the 90% confidence level. The case of pure standard oscillations, on the other hand, 
is ruled out at 90% confidence level. Note that only a fourfold rather than eightfold degeneracy is 
visible in fig. [4] because the transformation (13) only changes the sign of e™ and the octant of 623, 



which are not displayed here. 

The compatibility of the parameter region favored by MINOS with existing constraints depends 
strongly on whether atmospheric neutrino bounds are taken into account. Atmospheric neutrinos 
are very sensitive to non-standard matter effects since they can travel over very long distances 
inside the Earth before reaching the detector. Therefore, the MINOS-favored region of parameter 
space is in conflict with the limits derived in refs. |17t [T9] from MACRO and Super-Kamiokande 
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Figure 3: Comparison of MINOS data (black dots and error bars) to theoretical predictions including 
neutral-current NSI parameterized by the best fit point eq. (20) (blue dotted histograms) and charged 
current NSI parameterized by the best fit point eq. (21) (blue solid histograms). For comparison, the red 
dashed histograms show the theoretical prediction in the absence of neutrino oscillations, and the pink dash- 
dotted histograms represent the results of a two-flavor standard oscillation fit to the combined and 
data. 



data. While these bounds are numerically very strong, they have been derived in a two-flavor 
framework. It has been called into question whether they would still hold when standard and non- 
standard three-flavor effects are taken into account [161 131] . While we consider it highly unlikely 
that three-flavor effects would weaken the atmospheric neutrino bounds on NC NSI by the amount 
required to restore compatibility with the MINOS-favored parameter region, we cannot definitely 
rule out NC NSI as an explanation for the MINOS data at this time. 



3.2. Charged current NSI 



Let us now turn to the investigation of the charged current NSI introduced in sec. |2.2| We 
fit the MINOS data using a three-flavor analysis including the NSI parameters |e^J and arg(e^), 
but, as for the neutral current case, we find that three-flavor effects are small. In particular, the 



symmetries of the Vn survival probability found in sec. 2.2 for the two- flavor case are present also 



in the three- flavor framework. In particular, this means that there is a continuous family of best 



fit points (see eqs. (18) and (19), and the continuous generalization of (19) at the end of sec. 2.2). 
For reference, we here give one representative point from this family that is most consistent with 
existing bounds on 



0.12i = 0.12e 



0.5m- 



sm 



?23 



0.41 



Am 



32 



+2.74 x l(r d eV 2 



(21) 



with x 2 /dof = 15.9/15. The predicted MINOS event spectra at this point are shown as the solid 
blue histograms in fig. |3j 

In fig. [5, we show the allowed regions in the |e^|-arg(e^) plane determined by our fit. The 
continuous family of best fit points is indicated by the thin dash-dotted black curves. To the best of 
our knowledge, no bound on the vector-type v T -\i CC NSI considered here exists in the literature. 
(The bounds of order 0.1 derived in ref. [20J apply only to axial-vector operators.) However, the 
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Figure 4: Constraints on the parameter space of neutral current non-standard interactions in the /j,—t 
sector from MINOS and data. Each panel corresponds to a projection of the three-dimensional NSI 
parameter space along one of its axes. This means that in each panel we have marginalized over one of the 
NSI parameters (as well as the six standard oscillation parameters). The best fit points are indicated by the 
colored dots. The symmetries from eqs. (Ill and ( p~2[ ) are clearly visible in the plot, while the additional 
two-fold ambiguity eq. (13) is implicit. In the bottom panels, we explicitly indicate the parameter regions 
corresponding to a normal mass hierarchy (NH) and to an inverted mass hierarchy (IH), while in the top 
panel, NH and IH contours lie on top of each other. Exclusion limits from other experiments [TS] [20] are 
shown in gray. See text for caveats pertaining to bounds from atmospheric neutrino measurements. 



experimental limit on the branching ratio of the flavor-violating decay — > ^tt [32] can be 
translated into a bound of order 0.2 on |e^„| (and, in fact, also on the related coefficient |e^ r |) for 
vector and axial- vector type interactions (see appendix [B] for details and caveats). We expect that 
a bound could also be derived from lepton universality considerations in weak decays of parity-even 
hadrons, but that it would not be stronger than O(0. 1). Also, atmospheric neutrinos are sensitive 
to |e!jL|, but since CC NSI are not enhanced by the long baselines of atmospheric neutrinos, we 
expect these constraints to be relatively weak as well. 
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Figure 5: Constraints on the parameter space of charged current non-standard interactions between v T and 
muons from MINOS and data. We have used a full three-flavor fit, marginalizing over the standard 
oscillation parameters. The discrete and continuous symmetries from eqs. ( 18 1 and ( 19 1 are clearly visible in 
the plot. The thin dash-dotted black curves indicate the positions of the approximately degenerate best fit 
points. We explicitly indicate the parameter regions corresponding to a normal mass hierarchy (NH, blue) 
and to an inverted mass hierarchy (IH, red). See text for comments on existing constraints on CC NSI. 



4. TESTING THE NSI INTERPRETATION OF MINOS DATA IN FUTURE 

EXPERIMENTS 

To corroborate or refute the hypothesis of large non-standard interactions as an explanation 
for the apparent discrepancy between neutrino and anti-neutrino results in MINOS, it will be 
mandatory to gather more statistics in MINOS itself, and to look for possible NSI signals in future 
experiments like T2K and NOz^A. In the following, we will neglect neutral-current NSI since we have 
seen in sec. 3.1 that they are disfavored as an explanation for the MINOS data by atmospheric 
neutrinos. Instead, we will focus on CC NSI. We have computed the expected event spectrum 
in MINOS, T2K, and NOz^A assuming the values of the standard and non-standard oscillation 



parameters to be given by the MINOS best fit point eq. (21). We have then attempted a standard 
oscillation (no NSI) fit to this simulated data. If this fit is incompatible with the simulated data 
at a given confidence level, we say that the existence of a non-standard effect can be established 
experimentally at this confidence level. Our simulation of T2K follows [33l435| . while that of NO^A 
is based on [251 [36]. We include only the and disappearance channels. 

In fig. [6j we show the predicted discovery potential in MINOS, T2K, and NOz^A as a function of 
the integrated luminosity in neutrino mode and the integrated luminosity in anti-neutrino mode. 
We also indicate how the number of protons on target (pot) translates into a time of running at 
nominal luminosity (2.5 x 10 20 pot/year for MINOS, 6 x 10 20 pot/year for NOz^A, and 10 21 pot/year 
for T2K) . We see that optimal sensitivity is achieved if slightly more time is spent on running in anti- 
neutrino mode than on running in neutrino mode. This is easily understandable since the assumed 
NSI effect manifests itself mainly as an apparent discrepancy between neutrino and anti-neutrino 
results, while each data sample individually appears to be consistent with standard oscillations. 
Thus, to optimally probe the non-standard effect, the event numbers in the Vn and samples 
should not be too different. On the other hand, anti- neutrino cross sections are about a factor of 
3 smaller than neutrino cross section, so more time has to be devoted to v running to achieve this 
goal. Fig. [6] also shows that in order to improve the statistical significance of the anomalous effect 
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Figure 6: Discovery reach for charged current NSI parameters corresponding to the MINOS best fit point 



cq. (211 in MINOS (left), T2K (middle), and NOi/A (right). 



in MINOS itself, more anti- neutrino running is desirable since the experiment has already taken a 
lot of data in neutrino mode. By comparing the three panels of fig. [HJ we see that, as expected, 
the discovery potential of T2K is better than that of MINOS, while the best sensitivity is achieved 
in NOf A. After one year of nominal running in neutrino mode, NOi^A could confirm the existence 
of the non-standard effect at the 90% confidence level, while in anti-neutrino mode, even a few 
months would be sufficient to achieve that sensitivity. This can be understood by noting that, for 
the parameter values favored by MINOS, eq. (21), the two 0(|e^J) NSI terms in eq. (16) have 
opposite signs for neutrinos, but the same sign for anti-neutrinos. Therefore, the non-standard 
effect is stronger for anti- neutrinos. To achieve a 3a discovery in T2K or NO^A, neutrino and 
anti-running are required, with at least one year spent in each mode for T2K, or half a year for 
NOi/A. 

If the true values of the NSI parameters are different from the best fit point eq. (21), the 
discovery reach in future experiments can be altered significantly. This is illustrated in fig. [7j 
where we plot the \ 2 °f a standard oscillation fit to simulated data affected by CC NSI as a 
function of the running time at nominal luminosity. The widths of the colored bands correspond 
to the la uncertainty in the NSI parameters from fig. [5j Fig. [7] shows that if nature has chosen 
unfavorable NSI parameters, it will be very hard for T2K and NO^A to announce a discovery. On 
the other hand, for favorable parameter values a 3a effect could be detected after less than one 
year of nominal running even in T2K. 



5. DISCUSSION AND CONCLUSIONS 

We have seen that, in order to explain the tension between the Vu and z/„ event samples in 
MINOS using NSI, the NSI couplings would have to be rather large, almost of the same order 
as Standard Model weak interactions. While we have shown in sec. [3] that there are regions of 
parameter space still consistent with MINOS data and with constraints from other experiments, 
one should keep in mind that the effective operators generating the NSI should ultimately arise 
from an underlying renormalizable model. Model-dependent constraints, however, are usually much 
stronger than the model-independent bounds we have considered. 

For example, the most straightforward implementations of dimension 6 NSI operators, based 
on the introduction of new heavy tree-level mediator fields, are phenomenologically not viable 
because SU(2) invariance would dictate that large neutrino NSI realized that way would have 
to be accompanied by large non-standard effects in the charged lepton sector j37J EE]- There- 
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Figure 7: x 2 as a function of exposure time for a standard oscillation fit to (simulated) data affected by 
charged current NSI. The widths of the colored bands corresponds to varying the NSI parameters within 
the lcr allowed region preferred by MINOS (fig. [5]). We have assumed running at the indicated nominal 
luminosities, and we have assumed the running time to be equally divided into neutrino and anti-neutrino 
running. 



fore, such models are usually tightly constrained by rare decay searches [32]. NSI might arise 
from dimension 8 operators involving two Standard Model Higgs fields contracted with lepton 
doublets, so that after electroweak symmetry breaking, SU(2) breaking 4-fermion couplings can 
arise. However, dimension 8 operators of this type are typically accompanied by phenomenologi- 
cally problematic dimension 6 operators unless the coefficients of different operators obey certain 
cancellation conditions [38]. Further model-dependent constraints on neutrino NSI operators can 
come from electroweak precision tests such as muon g — 2 measurements, and from direct searches 
for possible mediators. All these constraints will typically force the mediators to be heavy (at least 
a few hundred GeV) or very weakly coupled. 

The latter possibility — neutrino NSI mediated by light (<C Mw), weakly coupled particles — is 
less well explored in the literature, so a scenario of this type could be responsible for the effects seen 
by MINOS. This is particularly interesting as models containing light new particles have recently 
received a lot of attention in the context of Dark Matter searches (see e.g. refs. [39— 44J). 

In conclusion, we have shown that the tension between the i/„ and disappearance data in 
MINOS — if it persists — could be explained by non-standard neutrino interactions. While neutral 
current interactions (non-standard matter effects) of the required magnitude are most likely ruled 
out by atmospheric neutrino constraints, a charged current operator leading to flavor- violating 
couplings between r-neutrinos and muons is not excluded. We have shown that such NSI can 
be tested in T2K and NOz^A, provided that the experiments are operated in neutrino and anti- 
neutrino mode. It remains an open question if and how NSI large enough to explain the MINOS 
results can arise from a renormalizable model. Along the way, we have derived the new constraints 
l e rJ < 0-20 and |e^ T | < 0.20 on flavor-violating vector or axial-vector type charged current non- 
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standard interactions in the /x— r sector by considering their loop-level contributions to the flavor- 
violating decay — > /i^vr . 
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Appendix A: GLoBES simulation of MINOS 

In this appendix, we provide details on the parameters of our MINOS simulation. We have used 
GLoBES [Ml [2S], with an implementation of NSI developed in refs. OSS], and with a MINOS 
experiment description based on [TJ I4"6l - I48| . The neutrino fluxes are taken from Monte Carlo 
simulations of the NuMI beam |49| . We use the same binning as the MINOS collaboration, but 
restrict our analysis to an energy window from 1-5 GeV for neutrinos, and from 1-8 GeV for anti- 
neutrinos. The reason for not including higher-energy neutrino data is that the fluxes available 
to us did not include the effect of unfocused high energy pions in the secondary beam. Detection 
efficiencies for Va, neutral current backgrounds, and the actual MINOS data were taken from pQ. 
The neutrino-nucleon scattering cross sections are based on the cross sections for water targets 
in |501 151] , and to account for the difference between the cross sections for water and those for 
iron, we adjust the overall normalization factor in our simulation in such a way that we optimally 
reproduce the predictions of the MINOS Monte Carlo simulation. We assume a Gaussian energy 
smearing function with a width oe given by Q.16E + 0.07y/E/GeV GeV for neutrino running, and 
by 0.155-E + O.llv^/GeV GeV for anti-neutrino running. These numbers were again optimized to 
reproduce the event rates predicted by the MINOS Monte Carlo simulation. The matter density 
along the neutrino trajectory is assumed to be 2.8 g/cm 3 . We have checked that the results of 
two-flavor standard oscillation fits to v„ and data agree well with the allowed regions obtained 
by the MINOS collaboration (see fig.[|]). 

In our fits, we use Gaussian prior terms to impose the la constraints sin 2 29\% < 0.1, sin 2 #12 = 
0.319 ± 0.023, and Am^ = (7.59 ± 0.30) x 10~ 5 eV 2 pi [53]. As systematic errors, we include 
independent 4% (3%) normalization uncertainties on the signal (background) rates for neutrinos, 
and 5% (5%) normalization uncertainties on the signal (background) rates for anti-neutrinos. 
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Figure 8: Allowed regions in the sin 2 2023 _ Am| 2 plane from separate two-flavor standard oscillation fits to 
MINOS and v^ L data. The thick contours and filled circles represent the 68% and 90% contours and best 
fit points from our fit, while thin lines and empty diamonds show the fit results obtained by the MINOS 
collaboration [T]. 
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Figure 9: (a) and (b) The loop diagrams used to constrain dj. and e^ T , respectively, (c) A similar diagram 
that does not have a logarithmic divergence. 



Appendix B: New loop bounds on non-standard interactions 



Here, we outline the calculation of the loop diagram fig. [9] (a) which we use to translate the 
experimental bound on the rare decay — > fi^ir into a constraint on |e!~J. We perform the 
calculation in unitary gauge and keep only the logarithmically divergent part. Indeed, it has been 
argued in refs. (18|, 120} [54"] that the finite and quadratically divergent terms cannot be calculated in 
a model-independent way, while the logarithmic divergence corresponds to the model-independent 
one-loop renormalization group running of the effective NSI operator from the UV completion 
scale down to Myy. (This argument is invalid if the UV complete theory leads to another effective 
operator having the same renormalization group running and exactly canceling the NSI operator. If 
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that case is considered, no model-independent bound on |e£„| can be derived.) Neglecting fermion 
masses, we find that the logarithmically divergent part of fig. [9] (a) is 

3V2G F e^,V ud a A r 

2 J{ l°g^[^][«7"ftu]. (BD 

where a is the electromagnetic fine structure constant, s w is the sine of the weak mixing angle, V u d 
is a CKM matrix element, Pl = (1 — 7 5 )/2, and A is the UV completion scale. Comparing eq. (Bl ) 



to the operator 2\f2GpV u d [t^PlVtWu^Pl^ responsible for the standard decay r 
using the QCD isospin symmetry between 7r° and -k^, and assuming logA/M^ ~ 0(1), we find 



d ^ / 2BR(r±-^/x±7r°) 4tt4 
^ - V BR(r±^7r±^ r ) 3a ' 1 J 

The factor 2 below the square root is due to the fact that the ir° can be interpreted as a state 
(uu — dd)/y/2 whose overlap with the final state uu from fig. § (a) is only l/y/2. With the 



experimental constraint BR(r =t — > [x^tt ) < 1.1 x 10 7 [32], eq. (B2) leads to the new bound 

l<J < 0-20. (B3) 
In full analogy to the above, we can also use fig. (b) to set a bound on e,'. 

I4r|<0.20. (B4) 

Note that diagrams like the one in fig. [9] (c) with dd in the final state do not contribute to these 
bound because they do not have logarithmic divergences. 



While we have derived eqs. (B3) and (|B4j) by assuming vector- type NSI, our calculations can 
be easily generalized to set equivalent bounds on axial- vector operators and on V — A operators, 
while no constraints can be derived for V + A type interactions since the diagrams in fig. [9] are 
suppressed by a neutrino mass insertion in this case. 



[1] P. Vahle (MINOS) (2010), Presentation at the XXIV International Conference on Neutrino Physics 
and Astrophysics (Neutrino 2010) in Athens, Greece; slides available at http://www.neutrino2010.gr/. 

[2] A. Dighe and S. Ray, Phys.Rev. D78, 036002 (2008), arXiv:0802.0121. 

[3] J. S. Diaz, A. Kostelecky, and M. Mewes (2009), 0908.1401. 

[4] G. Barcnboim and J. D. Lykken (2009), 0908.2993. 

[5] D. Choudhury, A. Datta, and A. Kundu (2010), 1007.2923. 

[6] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

[7] N. Engelhardt, A. E. Nelson, and J. R. Walsh (2010), 1002.4452. 

[8] Y. Grossman, Phys. Lett. B359, 141 (1995), hep-ph/9507344. 

[9] T. Ota, J. Sato, and N.-a. Yamashita, Phys. Rev. D65, 093015 (2002), hep-ph/0112329. 

[10] T. Ota and J. Sato, Phys. Lett. B545, 367 (2002), hep-ph/0202145. 

[11] N. Kitazawa, H. Sugiyama, and O. Yasuda (2006), hep-ph/0606013. 

[12] A. Friedland and C. Lunardini, Phys. Rev. D74, 033012 (2006), hep-ph/0606101. 

[13] M. Blennow, T. Ohlsson, and J. Skrotzki, Phys. Lett. B660, 522 (2008), hep-ph/0702059. 

[14] J. Kopp, M. Lindner, T. Ota, and J. Sato, Phys. Rev. D77, 013007 (2008), 0708.0152. 

[15] N. C. Ribeiro et al., Phys. Rev. D77, 073007 (2008), 0712.4314. 

[16] M. Blennow, D. Meloni, T. Ohlsson, F. Terranova, and M. Westerberg, Eur. Phys. J. C56, 529 (2008), 
0804.2744. 

[17] N. Fornengo, M. Maltoni, R. T. Bayo, and J. W. F. Valle, Phys. Rev. D65, 013010 (2001), hep- 
ph/0108043. 



16 



[18] S. Davidson, C. Pena-Garay, N. Rius, and A. Santamaria, JHEP 03, Oil (2003), hep-ph/0302093. 

[19] M. C. Gonzalez-Garcia and M. Maltoni, Phys. Rept. 460, 1 (2008), 0704.1800. 

[20] C. Biggio, M. Blennow, and E. Fernandez-Martinez (2009), 0907.0097. 

[21] W. A. Mann, D. Cherdack, W. Musial, and T. Kafka (2010), 1006.5720. 

[22] E. Akhmedov and T. Schwetz (2010), 1007.4171. 

[23] J. Heeck and W. Rodejohann (2010), 1007.2655. 



[24] K. Nishikawa et al. (T2K) (2006), URL |http : // j -pare . jp /NuclPar t/pac_0 606/pdf / 

pll-Nishikawa.pdf . 

[25] D. S. Ayres et al. (NOz^A) (2004), hep-ex/0503053. 

[26] R. Gandhi, P. Ghoshal, S. Goswami, P. Mehta, and S. Uma Sankar, Phys. Rev. Lett. 94, 051801 (2005), 
hep-ph/0408361. 

[27] P. Astier et al. (NOMAD), Nucl. Phys. B611, 3 (2001), hep-ex/0106102. 

[28] P. Huber, M. Lindner, and W. Winter, Comput. Phys. Commun. 167, 195 (2005), hep-ph/0407333, 

URL http : // www . mpi-hd . mpg . de/ -globes 

[29] P. Huber, J. Kopp, M. Lindner, M. Rolinec, and W. Winter, Co mput. Phys. Commun. 177, 432 (2007), 

hep-ph/0701187, URL http://www.mpi-hd.mpg.de/~globes 

[30] M. Maltoni and T. Schwetz, Phys. Rev. D68, 033020 (2003), hep-ph/0304176. 

[31] A. Friedland and C. Lunardini, Phys. Rev. D72, 053009 (2005), hep-ph/0506143. 

[32] C. Amsler et al. (Particle Data Group), Phys. Lett. B667, 1 (2008). 

[33] P. Huber, M. Lindner, and W. Winter, Nucl. Phys. B645, 3 (2002), hep-ph/0204352. 

[34] Y. Itow et al. (2001), hep-ex/0106019. 

[35] M. Ishitsuka, T. Kajita, H. Minakata, and H. Nunokawa, Phys. Rev. D72, 033003 (2005), hep- 
ph/0504026. 

[36] T. Yang and S. Woijcicki (NOi/A) (2004), Off-Axis-Notc-SIM-30. 

[37] S. Antusch, J. P. Baumann, and E. Fernandez-Martinez (2008), 0807.1003. 

[38] M. B. Gavela, D. Hernandez, T. Ota, and W. Winter (2008), 0809.3451. 

[39] A. L. Fitzpatrick, D. Hooper, and K. M. Zurek (2010), 1003.0014. 

[40] D. Hooper, J. Collar, J. Hall, and D. McKinsey (2010), * Temporary entry *, arXiv:1007.1005. 

[41] C. Arina and M. H. G. Tytgat (2010), 1007.2765. 

[42] J. Lavalle (2010), 1007.5253. 

[43] D. T. Cumberbatch, J. A. Guzik, J. Silk, L. S. Watson, and S. M. West (2010), 1005.5102. 

[44] E. Kuflik, A. Pierce, and K. M. Zurek (2010), 1003.0682. 

[45] J. Kopp, M. Lindner, and T. Ota, Phys. Rev. D76, 013001 (2007), hep-ph/0702269. 

[46] E. Abies et al. (MINOS) (1995), FERMILAB-PROPOSAL-0875. 

[47] P. Huber, M. Lindner, M. Rolinec, T. Schwetz, and W. Winter (2004), hep-ph/0403068. 

[48] B. V. M. Diwan, M. Messier and L. Wai, NUMI-L-714 (2001). 

[49] M. Bishai (2010), private communication. 

[50] M. D. Messier (1999), UMI-99-23965. 

[51] E. A. Paschos and J. Y. Yu, Phys. Rev. D65, 033002 (2002), hep-ph/0107261. 

[52] M. Maltoni and T. Schwetz (2008), 0812.3161. 

[53] M. C. Gonzalez-Garcia, M. Maltoni, and J. Salvado (2010), 1001.4524. 

[54] C. Biggio, M. Blennow, and E. Fernandez-Martinez, JHEP 03, 139 (2009), 0902.0607. 



